Dark matter-wave solitons in the dimensionality crossover 
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We consider the statics and dynamics of dark matter-wave solitons in the dimensionality crossover 
regime from 3D to ID. There, using the nonpolynomial Schrodinger mean-field model, we find that 
the anomalous mode of the Bogoliubov spectrum has an eigenfrequency which coincides with the 
soliton oscillation frequency obtained by the 3D Gross-Pitaevskii model. We show that substan- 
tial deviations (of order of 10% or more) from the characteristic frequency u z /\/2 (cj 2 being the 
longitudinal trap frequency) are possible even in the purely ID regime. 



The experimental realization of lower-dimensional 
Bose-Einstein condensates (BECs) in highly anisotropic 
traps [3, [|J has inspired many studies devoted to the be- 
havior of such systems in the dimensionality crossover, 
i.e., from three-dimensions (3D) to one-dimension (ID). 
Particularly, it has been shown that fundamental prop- 
erties of BECs, such as the chemical potential, speed of 
sound, collective oscillations, change significantly as the 
dimensionality is reduced from 3D to ID 0, 0|- These 
regimes, as well as the crossover between them, can be 
described [4| by the dimensionless parameter N£la/a±, 
where N is the number of atoms, Q = lu z /lu±_ is the ra- 
tio of longitudinal and transverse trapping frequencies, 
a is the scattering length, and a± = y / nJrnuT[ (m is 
the atomic mass). In particular, if the dimensionality 
parameter is Nila/a± ^> 1, the BEC locally retains its 
original 3D character and its ground state can be de- 
scribed by the Thomas-Fermi approximation in all di- 
rections. On the other hand, if Nfla/a± <§; 1, excited 
states along the transverse direction are not energetically 
accessible and the BEC is effectively ID; in such a case, 
the transverse wave function is described by the ground 
state of the radial harmonic oscillator, whereas the lon- 
gitudinal wave function is governed by an effectively ID 
Gross-Pitaevskii equation (GPE) Q. Of particular rele- 
vance are effectively ID mean- field models [fall! that have 
been developed to describe the axial dynamics of "cigar- 
shaped" BECs in the ID (and 3D) regimes, as well as in 
the dimensionality crossover. Importantly, there exist re- 
cent experimental results [8|| which have been accurately 
described by the model of |7|. 

It is then natural to expect that, apart from the 
ground-state properties of BECs, the dimensionality 
should significantly affect the stability and dynamical 
properties of the excited BEC states as well, such as 
the dark solitons in repulsive BECs. The first exper- 
iments reporting the observation of dark matter-wave 
solitons were performed with quasi-spherical 0, [l(| or 
cigar-shaped [lj| traps. In all cases, dark solitons were 
prone to instabilities, such as non-uniformity induced 
dynamical instability (leading to a U-shaped deforma- 
tion of the soliton propagation front) thermal in- 



stability (decay due to the interaction with the thermal 
cloud) [111 ], or snaking instability (leading to soliton de- 
cay into vortex rings) • From the theoretical point of 
view, the former two instabilities were analyzed in [T3 |. 
while the snaking instability (which occurs in the higher- 
dimensionality setting [3]) was analyzed in [iH Il5l. flB|. 
In the latter works, a detailed study of the Bogoliubov- 
de Gennes (BdG) equations in 2D and 3D revealed the 
emergence of complex eigenvalues in the excitation spec- 
trum and their connection to oscillatory dynamical in- 
stabilities (including the sna king instability). Moreover, 
it was demonstrated [3, 15 , Hit [13] that there exists a 
so-called "anomalous mode" in the excitation spectrum, 
which has negative energy. The existence of the anoma- 
lous mode indicates that dark solitons (and vortices [lil ]) 
are thermodynamically unstable and, in the presence of 
dissipation, the system is driven towards configurations 
with lower energy; this scenario is also often referred to 
as energetic instability (l~9j . Importantly, the anomalous 
mode frequency is directly related to the dark soliton 
oscillation frequency uj osc = u> z / y/2 [3 E3] ■ This oscil- 
lation frequency has been obtained upon analyz ing the 
ID GPE using different analytical approaches |2fJ for 
f2 -C 1, and assuming the Thomas- Fermi (TF) approx- 
imation for the longitudinal BEC background on which 
the dark soliton exists. However, it is important to stress 
that the above result is not particularly relevant to either 
the first dark soliton experiments [^, [T(| [ll[ , or to the re- 
cent experimental studies of [2l|, because of the fact that, 
in all cases, the condensates were actually in the purely 
3D regime. 

In this work, we revisit this problem and study the 
statics and dynamics of dark solitons in BECs confined 
in highly elongated traps (with Q <C 1) in the dimen- 
sionality crossover regime, i.e., for Ntta/a± ~ 1; our 
analysis also applies in the case of quasi-lD small BECs 
in which the TF approximation for the axial direction 
is not valid. We thus aim to provide a complete pic- 
ture bridging the one- and multi-dimensional dynamics 
of dark solitons. We will show, in particular, that in the 
crossover regime the anomalous mode frequency resulting 
from the BdG analysis of the nonpolynomial Schrodinger 
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equation (NPSE) Q estimates accurately the soliton os- 
cillation frequency resulting from the 3D GPE. Moreover, 
we will demonstrate that, in this regime, the frequencies 
of the anomalous mode and the soliton oscillation coin- 
cide but may significantly differ from the value lu z /\/2 
(deviations may easily exceed 10%). 

Let us start our analysis considering the example of a 
87 Rb condensate containing N ~ 2100 atoms and con- 
fined in a trap with frequencies wj_ = 20u> z = 2tt x 200 
Hz (here ft = 0.05). In this case, the dimensionless pa- 
rameter NQa/a± = 0.815, a value in the dimensionality 
crossover region, and thus the ID GPE is not valid (the 
transverse wave function deviates from the ground state 
of the transverse harmonic oscillator). Thus, to calculate 
the oscillation frequency of a dark soliton in such a set- 
ting, one needs to consider the fully 3D GPE, which can 
be expressed in the following dimensionless form: 
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where the density \ip\ 2 , length, time and energy are re- 
spectively measured in units of Airaa 2 ^, a±, ui7 and thu±. 
We use as an initial condition a dark soliton initially 
placed at z = 2 on top of the ground state of the sys- 
tem, which is found by a relaxation technique (i.e., us- 
ing imaginary time integration to find the ground state, 
as well as a quiescent dark soliton at the trap center, 
and then displacing the soliton at the desired location 
to set the soliton into motion). Then, for a normal- 
ized 3D chemical potential /i^d — 1.625, the following 
result is obtained (see Fig. 1): The soliton oscillates 
with a frequency uj osc = 0.746ft, although the system is 
"highly anisotropic" and one would expect an oscillation 
frequency of ft/\/2 = 0.7070 (in this case the deviation 
from the prediction of [20j is of order of » 6%). 

Instead of using the fully 3D GPE model, it is worth 
investigating if the above value of the soliton oscillation 
frequency can be obtained in the framework of a ID mean 
field model. As the ID GPE is not applicable in this 
example, we consider another relevant candidate, namely 
the NPSE of Ref. 0]. The latter, can be expressed in 
the following dimensionless form, 



dt 



\cP_ 



-tfz 2 



3|# 



2(1 + |0| 2 )V2 



(2) 



where <fi is the normalized longitudinal wavcfunction, and 
the density |</>| 2 is measured in units of 2aN (length, time 
and energy are measured in the same units as in the 3D 
GPE). Note that this model is reduced to the ID GPE 
in the low density limit of \<fi\ 2 <C 1 (in our units), which 
is not relevant to the case example of N£la/a± = 0.815. 
We have numerically integrated Eq. ([2]) with an initial 
condition corresponding to longitudinal density profile of 
the 3D GPE case and have found that the soliton oscilla- 
tion frequency is uj osc = 0.744ft, approximately the same 
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FIG. 1: (Color online) (a) Contour plot of the initial density 
(in the x = plane) of a BEC, confined in an harmonic trap 
with fl = 0.05, and with a dark soliton placed at z = 2. The 
dimensionality parameter is NQ.ot/a±_ = 0.815. (b) Spatio- 
temporal evolution of the above density along a cross-section 
at r = 0. The dark soliton oscillates with frequency cu osc = 
0.746ft > Q/y/2. The dotted line across the soliton trajectory 
corresponds to the prediction obtained by the NPSE model. 
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FIG. 2: (Color online) (a) The six lowest characteristic eigen- 
frequencies of the Bogoliubov excitation spectrum: Two of 
them, located at the origin, correspond to the Goldstone mode 
due to the phase invariance, the ones at ±ft = ±0.05 to the 
Kohn mode, while the indermediate ones, at ±0.74457, to the 
anomalous mode. The dashed lines indicate the location of 
the ±J7/v2 prediction, (b) The eigenfunctions u (solid line) 
and v (dotted line) of the anomalous mode. 



with the one obtained by the 3D GPE (w osc = 0.746ft). 
This result, depicted by the dotted line in Fig. [IJb), 
indicates that the NPSE may accurately predict the soli- 
ton oscillation frequency in the dimensionality crossover 
regime. 

Following the above argument, we now proceed with 
the investigation of the excitation spectrum of the 
NPSE's dark soliton state. This way, we first find a real, 
stationary soliton state 4>ds(z) using a fixed point algo- 
rithm (a standard Newton- Raphson method [22j). Then, 
considering small perturbations of this state of the form, 
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we obtain from Eq. ([2]) the following BdG equations 
(valid to leading order in the small parameter e): 

">« = + /(&)]« + fl($>, (4) 

-uv = [#-/* + /($.)]*> + fl($>, (5) 



where H = -(l/2)S| + (l/2)f2 
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and g{4>\ s ) = 4 ^ d3 ^9 ^3/2 • The above equations provide 

the eigenfrequencies w = w r + iuJi and the amplitudes 
u and t> of the normal modes of the system. Note that 
due to the Hamiltonian nature of the system, the eigen- 
frequencies of the Bogoliubov analysis appear in pairs 
(or in quartets, if complex); thus, the solution of BdG 
equations with frequency uj represent the same physical 
oscillation with the solution with frequency — oj [231 ]. 

Among the various eigenfrequencies, we focus on the 
three smallest magnitude pairs in Fig. |2[a): One of them 
is at the origin of the spectral plane (u} r ,u>i) reflecting 
the phase invariance of the NPSE model. The respective 
eigenfunction is the so-called Goldstone mode and does 
not result in any physical excitation (oscillation) of the 
system. On the other hand, the solutions with eigenfre- 
quencies u> = fl = 0.05 correspond to the so-called dipole 
(or Kohn) mode, representing the motion of the center of 
mass of the system which oscillates with the frequency of 
the harmonic trap. Finally, of particular interest are the 
solutions with eigenfrequencies u — 0.744S1 correspond- 
ing to the so-called anomalous mode. For the latter, the 
integral of the norm x energy product, J u>(\u\ 2 — \v\ 2 )dz 
(in our units), is negative [23(, indicating the energetic 
instability of the dark soliton discussed above. Impor- 
tantly, this eigenfrequency coincides with the oscillation 
frequency of the dark soliton as obtained by the NPSE 
model-and in accordance with the result of the 3D GPE. 
The eigenfunctions u and v of the anomalous mode are 
shown in Fig. El^b). As seen, they are localized within 
the notch of the dark soliton [l7| . 

We now aim to investigate if the above agreement 
is generic, i.e., whether the BdG analysis of the NPSE 
model indeed provides an accurate estimation of the soli- 
ton oscillation frequency in the dimensionality crossover 
regime. In this connection, we fix the normalized har- 
monic trap strength to its previous value (CI = 0.05), and 
find the soliton oscillation frequency varying the param- 
eter Na/a±. The results are presented in Fig. [3] which 
summarizes the main findings of the present work. The 
obtained normalized frequency u/Cl is shown: the solid 
line depicts the anomalous mode frequency obtained by 
the BdG analysis of the NPSE model, while the stars de- 
pict the soliton oscillation frequency obtained by a direct 
numerical integration of the fully 3D GPE. The differ- 
ent dimensionality regimes, are also shown in Fig. [3l 
and are defined as follows: The ID regime corresponds 
to fi3D < hw±_ or < 1 in our units, i.e., to 

Na/a± < 2 (or Nfla/a± < 0.1) as per our normal- 
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FIG. 3: (Color online) Normalized frequency lu/Q. as a func- 
tion of Na/a±: Solid (dashed) line corresponds to the anoma- 
lous mode frequency obtained by the BdG analysis of the 
NPSE (ID GPE), and stars denote the oscillation frequency 
of the dark soliton obtained by the 3D GPE. The vertical 
dotted lines depict the different dimensionality regimes. 



izations; on the other hand, the dimensionality crossover 
regime corresponds to the region around Na/aj_ ~ 20 
(or Nfla/a± as 1), while the 3D regime corresponds to 
the limit Na/a ± > 20 (or NCla/a x > 1). 

In the ID regime, the NPSE model is reduced to the 
ID GPE and, as a result, the anomalous mode frequency 
obtained in the framework of the NPSE coincides with 
the one obtained by the BdG analysis of the ID GPE (see 
dashed line in Fig. [3]). This agreement ceases to exist for 
Na/a± > 2, which is a clear indication that the system 
enters the dimensionality crossover regime. Therefore, 
the asymptotic value of fl/y/2 depicted by the dotted line 
in Fig. [3]for Na/a± 3> 2 is quantitatively irrelevant: The 
correct result for the oscillation frequency in the crossover 
regime is provided by the 3D GPE (stars) and is accu- 
rately approximated by the NPSE (solid line). Notice 
that in the limit Na/a± — > (which corresponds to the 
linear Schrodinger equation), one obtains lu/Q — > 1, i.e., 
the anomalous and Kohn mode frequencies coincide, in 
accordance to the prediction of Ref. 
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In the dimensionality crossover regime, and particu- 
larly in the interval 2 < Na/a± < 20, the anomalous 
mode frequency obtained by the BdG analysis of the 
NPSE almost coincides with the oscillation frequency of 
the dark soliton obtained by the 3D GPE. However, as 
shown in Fig. [31 as the system approaches the 3D regime 
(Na/a± 3> 20) the NPSE underestimates the frequency 
obtained by the 3D GPE. This is a consequence of the 
fact that the NPSE uses a gaussian ansatz to describe 
the transverse wavefunction, rather than the TF profile 
which is relevant to the 3D regime. On the other hand, 
it is important to note that for (i 3 d > 2.4 14, [lj|, corre- 
sponding to Na/a± as 60, the system becomes dynami- 
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cally unstable through the emergence of complex eigen- 
frequencies in the excitation spectrum of the 3D GPE, 
related to the onset of oscillatory and snaking instabili- 
ties of dark solitons 0, [Til, [l6| . Although these types 
of instabilities were indeed observed in our 3D simula- 
tions (results not shown here), the BdG equations of the 
NPSE did not predict such complex eigenvalues. Thus, 
the above results suggest a note of caution: the BdG 
analysis of the NPSE model can indeed predict accurately 
the oscillation frequency of dark solitons in the dimen- 
sionality crossover regime, but, given its freezing of the 
transverse directions into their ground state, it cannot 
capture the emergence of the pertinent instabilities oc- 
curing in the fully 3D regime. 

We have also performed the BdG analysis of both the 
NPSE and ID GPE models for other values of the nor- 
malized trap frequency £1. As seen in Fig. [4l where 
cu/fl is shown as a function of the dimensionality param- 
eter NQa/a±, larger (smaller) values of £1 yield larger 
(smaller) anomalous mode frequencies. It is worth men- 
tioning that Fig. 2] is merely devoted to the ID regime 
(the dimensionality parameter does not exceed the value 
0.3). There, as the trap frequency is decreased (towards 
the limit of £1 -C 1 considered in Refs. [1(3]), the anoma- 
lous mode frequency approaches the value of Q/v2 de- 
picted by the dashed line in Fig. 0J Thus, in the ID 
regime, and for £1 <§; 1, the BdG analysis of the NPSE re- 
covers the soliton oscillation frequency obtained asymp- 
totically in the framework of the ID GPE [201 ] . 

An important observation stemming from our analy- 
sis is that even in the ID regime, the soliton oscillation 
frequency may have a substantial difference from £1 / y/2. 
For example, let us consider two different cases (both 
for a 87 Rb BEC), one with number of atoms N ~ 5000 
and trap frequencies lu± — 200cj 2 = 2ir x 200Hz, and 
one with N ~ 1000 and ui±_ = 10lu z = 2tt x 70Hz. In 
the former case, our analysis predicts that the soliton os- 
cillation frequency will be uj osc = 0.718£1, differing only 
1.5% from the value of fl/y/2, while in the second case 
uj osc = 0.772£1, differing 10% from £1/V2. Such strong 
deviations from the asymptotic prediction should be di- 
rectly accessible to current experimental settings. 

In conclusion, we have studied dark-matter wave soli- 
tons in Bose-Einstein condensates and how their dynam- 
ics is affected by dimensionality in settings of experimen- 
tal interest/accessibility. We have first shown that in a 
highly anisotropic system at the crossover between 3D 
and ID behavior, the dark-soliton oscillation frequency 
is larger than the usually quoted value of fl/y/2. The 
deviations are predicted by the ID mean-field NPSE 
model, whose excitation spectrum reveals an anomalous 
mode eigenfrcquency identical to the soliton oscillation 
frequency in the dimensionality crossover regime. Limi- 
tations of the NPSE mdel, concerning its validity towards 
the fully 3D regime (where dynamical instabilities of dark 
solitons may also come into play) were also discussed. 
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FIG. 4: (Color online) Normalized frequency uj/Q. of the 
anomalous mode obtained by BdG analysis of the NPSE 
(solid lines) and the ID GPE (dotted lines), as a function of 
Ntla/a± for £1 = 0.1,0.05,0.005. The dashed line indicates 
the fi/\/2 anomalous mode eigenfrequency. 



Importantly, our analysis demonstrates that, even in the 
purely ID regime, deviations of the soliton oscillation 
frequency from the standardly used theoretical value of 
of order of 10% (or even more) are possible and 
should be observable in current experimental setups. 
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